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We study a (2+l)-dimensional system that can be viewed as an infinite number of 0(3) cr-fields 
coupled by a nearest-neighbour Heisenberg-like interaction. We reduce the field equations of this 
model to an integrable system that is closely related to the two-dimensional relativistic Toda chain 
and the Ablowitz-Ladik equations. Using this reduction we obtain the dark-soliton solutions of our 
model. 
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' 1. Introduction. 



The model considered in this paper can be viewed as a generalization of the classical 0(3) 
(T-model in two dimensions, described by the Hamiltonian function 



£ = £[(r] = j dxdy (Ver,Ver) (1.1) 



where a is a three- component vector of unit length, 

(*,<r) = l. (1.2) 
and braces denote the standard scalar product. The energy of our system is given by 

n = J2 £ [°n]+U in t (1.3) 

n 

with nearest-neighbour interaction 

Hint=\Y. E U -P (1-4) 

n p= n ±l 



1 
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of the Heisenberg type: 



U np =U[a n ,a p ] = J dxdy F f(cr n ,cr p )\ . (1.5) 

The models of this type can appear, for example, in the studies of the lamellar (graphite- 
like) magnetics when the spin interaction inside one layer can be described in the framework 
of the Landau-Lifshitz theory with effective Heisenberg interaction between adjacent layers. 
The stationary structures of our system are governed by the (2+l)-dimensional equation 

0, (<5<x n ,cr n ) = 0. (1.6) 



5<j t 

In what follows we use a function F which is peculiar to integrable nonlinear mathematics 
(see e.g. [HE]), 

F(x) = 5 2 ln(l + x). (1.7) 
The resulting equations are given by 

[Aer n ,er n ] = — ^ fnp[c p ,(T n ] (I- 8 ) 

p=ra±l 

where 

2 

1 + [cr n , cr p ) 

The factor g 2 can be eliminated by rescaling the coordinates, so we take 

5 = 4 (1.10) 
and write the central equation of our study as 

\ [Aa n ,a n ] = ^ fnp[o-p,CTn\- (1-11) 

p=n±l 

In the following sections, after re-parametrization of ([l.lip . we split it in Sec. [2] into a 
first-order system, bilinearize it (Sec. [3]) and derive the dark-soliton solutions (Sec. |3J). 

2. Parametrization and splitting 

Using the vector-matrix correspondence 



where cr J (j = 1, 2, 3) are the Pauli matrices 



io/' u or vo — i 
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and introducing complex variables 

z = x + iy, z = x — iy (2.3) 

one can rewrite Eq. (jl.lip as 

§n,§n] = Yl fnpi^n] (2.4) 
p=n±l 



with d = d/dz, d = d/dz and 



fnp ~ 1 + ltrSnSp- (2 ' 5) 



In what follows we use the parametrization of the vectors cr n based on the presentation of 
the matrices § n in the form 

§ n = ^"V 3 ^ n . (2.6) 

Using the invariance of this representation with respect to transformations \& n — > B n \I / n 
with arbitrary diagonal matrices D n one can choose 



which leads to 



g _ 1 /l + B n C n 2B n , _ 
1 — B n C n \ —2C n —1 — B n C n 



Calculating ddS n and /, 



j = (1 - B n C n ) (1 - B p Cp) ^ ^ 

p (1 — B n C p ) (1 — C n B p ) ' 



one comes to the following system of equations: 

f A n Cn = Y n (B n+ i — B n ) — l" n _i (B n - B n - 

I A n £n = (Cn+1 — C n ) — Yn— 1 (Cn — C n _ 



where 



and 



(2.10) 



= 1 „ ^ , (2-11) 



£^ = c*d5 n + 2^ n (55 n ) (&B n ) C n , (2.12) 
= d<9C n + 2A n B n (dC n ) (BC n ) (2.13) 
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The crucial step of our proceeding is the following ansatz: we split the above system 
into two first-order ones, 



and 



iA n dB n = Z n -\ (B n - B n -i) 
iA n dC n = Z n (C n +i — C n ) 



-iA n d_B n = Z n (B n+l -B n ) 
-iA n dC n = Z n —i (C n — C n _i) . 



(2.15) 



(2.16) 



By direct calculations one can show that this can be done provided we can find the func- 
tions Z n and Z n that (i) make (|2.15p and (|2.16p compatible and (ii) lead to (|2.10p . It is 
demonstrated in the appendix that the functions Z n and Z n that meet these conditions can 
be chosen as 

Z n = (Y n , Z n = (Y n (2.17) 
where £ and £ are arbitrary constants related by 

CC = 1 (2.18) 
To summarize, one can obtain a large number of solutions of (j2.4j) by solving the system 

n - (B n - B n -{) 

(2.19) 



idB n = C 
idC n = C 



1 — B n -\C n 

• ~ B n C n 
1 — B n C n +\ 



(B n - B n -i) 
{Cn+l ~ Cn) 



and 



-idB n = C / * nCn r (B n+1 - B n ) 

-iBc n = c : 1 ~ B B t r n (c n - c n _x) 

1 — i( n G n _i 



(2.20) 



Before proceed further, we would like to give some comments on this system. After 
introducing new variables, 



B n — 1/ C n , 

Eqs. (|2.19|) . (|2.20|) can be cast into the Hamiltonian form 

2 dH 



(2.21) 



idB r , 



B n — B r 



idB n = [B n — B n 



dB n 

2 dH 



(2.22) 



dB n 



with 



h = c jr in Bn ~ Bn 

n=-oo B n - B n+ i 



(2.23) 
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and 



idB n 
—idB r , 



B n — B n 



B n — B n 



2 8H 

2 oh 



(2.24) 



with 



H = C In 



dB n 

B n — B n 
B n — B n _i 



(2.25) 



and can be identified with the (Xi,Yi) equations (with a(u,v) = (u — v) 2 ) from the list of 
the paper by Adler and Shabat [3]. 

At the same time both B n and C n solve the (2+l)-dimensional version of the 
Ruijsenaars-Toda lattice [HE] 



dBU n + (dU n ) [dU n ) 



u, 



n+1 



u„ 



Un ~ U, 



n—1 



0. 



(2.26) 



Note that equations (|2.26p are different from (and complementary to) the Ruijsenaars-Toda 
lattice (i?i) that appears in a natural way in the framework of [3]. 

Finally, calculating from Eqs. (|2.15p . (|2.16p derivatives of the functions f n defined by 



fn — /n,n+l 



(1 — B n C n ) (1 — B n+ iC n+ i) 



(1 — B n C n+ \) (1 — B n+ \C n ) 
one can demonstrate that these functions satisfy 

dd In /„ = f n+1 - 2f n + / n _i. 



(2.27) 



(2.28) 



Thus one can see the relationship of the model discussed in this paper with the famous 
two-dimensional Toda lattice. 



3. Bilinearization. 

To bilinearize Eqs. (|2.19p . (|2.20p we introduce p n , f n , f n and a n by 

Pn-l ~ °r, 



B n 



Cm 



T n -1 



and another set of tau-functions by 



and 



iDpn-l ■ T n = a Pn-lTn 

iD f n _i ■ a n = a r n -\a n 



-iDPn-l -T n = a T n _ip n 

-iDf n -\ -& n = a (J n ^\T n 



(3.1) 



(3.2) 



(3.3) 



where a and a are constants, D and D are the Hirota's bilinear differential operators, 
Du ■ v = (du) v — u (dv) and Du ■ v = (du) v — u (dv) . Now, to finish the bilinearization 
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of our equations, we impose the restrictions 

f PnTn - Pn-lTn+1 = P p n T n 

I f n a n - f n _i<7 n+ i = (3 T n (j n 



(3.4) 



and 



f T n _lf n + Pn-l&n = 1 A T n -lT n 
f n f n + Pn-lVn+1 = 1 B T% 

_ f n _if n+ i + p n a n = 7 C tI 



(3.5) 



where j3, 7 , ^ B and ^ c are again some constants. It can be shown by direct calculations 
that Eqs. (|5^]> -([53 ]l imply that B n , C n satisfy Eqs. (I2TT9]) . (I2T20]) . Indeed, noting that Eqs. 
(|3.4p and (|3.5p are nothing but 



B n +i — B n = (3 



Cn+l ~ C n — (3 — 



PnTn 
T n T n +l 
Tn&n 



and 



1 - B n C n = 7 
1 — B n C n+ \ = 7 
1 — B n+ iC n = 7 ( 



Tn—l^n 



A Tn—lTn 
_2 



TnTn 



Tn-lTn+l 



calculating f n , 



(7 ) T n _iT n+ l 



7 7 



and substituting the above formulae into f)3. 2[) and ()3.3[) one can obtain 

1 — B n C n 



idB n = T 1 



iOC n = r 



c 



1 — B n -iC r 

1 B n C n 



(B n — B n _\ 



1 — B„.C, 



(Cn+l ~~ C n ) 



and 



ias n = r B 1 gwC " (s n+1 - s n ) 



1 — B n C n -\ 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



(3.11) 

(3.12) 
(3.13) 

(3.14) 
(3.15) 
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where 



7 B a 



-.c 



7 a 



c _ 

7 a 



f 



7 a 



7 A /3' 7 A /T 7^/3' 7 A /3' 

Thus, to finish solution of our problem one has impose the condition 



c, 



c = vc 



(3.16) 



(3.17) 



In this way we have reduced Eqs. (|2.19[) . (|2.20p . and hence Eqs. (|2.4p . to the set of 
the bilinear equations f|3 . 2 p — f|3 . 5 [) . An important question that arises now is the question 
about compatibility of this system. We do not present here an explicit proof of the fact that 
Eqs. (|3.2|) - (|3.5p are compatible because (i) we present (in the next section) their explicit 
solutions and (ii) show their relation to a well-known nonlinear compatible system - the 
Ablowitz-Ladik hierarchy (ALH) [6j. To do the latter let us consider the matrix 



1 



Tn-l 



T n Pn-1 

-<7 n f n _l 



Calculating its determinant, 



and inverse one can obtain 



with 



det $ n = Y 



Tn-l 



and 



n+1 



Un = — 









Tn 






1~n 





(3.18) 



(3.19) 



(3.20) 



(3.21) 



with 



and 



id® n = V n <5> n 



V n 



idhi 



Tn 



a Pn-l&r, 



+ 

T n -1 7^ T n ^iT n 
1 A T n 



id$ n = V n <S> n 



a Pn-1 
7 A Tn-l 



id In Tn 1 



T n -1 



7 



0_ Pn-\On 
1 f n -lT n / 



T, 



Vn 



d In r 



Tn-l 

a a n -\ 



1 A T n f n 



a Pn 
7 A T n 



1 A Tn-l 



.s, T n -\ a p n -\a n -\ 
i d m 1 -r — 



Tn-l 



7^ T n _iT„ 



1 / 



(3.22) 



(3.23) 



(3.24) 
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Inspecting (|3,20p - (|3.24p one can conclude, after eliminating the unnecessary constants, 
introducing 

q n = — , r n = — (3.25) 

and making some simple gauge transformations, that (|3.20p with (|3.2ip is nothing but the 
spectral problem of the ALH whereas Eqs. (|3.22p with (|3.23p and (|3.24p describe its first 
positive and negative flows. So, bilinear equations (|3.2p - (|3.5p belong to the ALH. This 
leads to two important results: (i) they are compatible and (ii) we can use already known 
solutions for the ALH to get solutions of our equations. 

To expose the inner structure of Eqs. (|3.2|) - (|3.5p and to make the following formulae 
more readable it seems useful to introduce instead of the triplet p n , r ra and a n an infinite 
set of tau- functions r™, 

Pn = T~ 1 , T n = r®, (y n = r\. (3.26) 
In new terms equations (|3.2p - (|3.5p become 



LUT n-\ T n 


„ ^fn—X^rn 
= « T n~l T n 


(3.27) 


—if)f m ~ X ■ T m 
LUT n-\ T n 


= OT m ,T m_1 


(3.28) 


-m— l~m ~m—\±.m 
'n 'n n—l n+1 


= T m ~ l T m 


(3.29) 



for m = 0, 1 and 



C-iC + C-tc +1 = r C-i^ (3-30) 

CC + C-l 1 ^ 1 = 1 B (C) 2 (3.31) 
+ -™_ iC+i = 7 c (C)2 (3 _ 32) 

for m = 0. These equations are a part of the generalized ALH |7] and can be solved without 
imposing restrictions on to, for — oo < to < oo. 



4. Dark solitons. 

4.1. Dark solitons of the ALH. 

Here we would like to present some basic formulae describing the dark-soliton solutions of 
the ALH that we then use to obtain solutions of our problem. 

The dark solitons for the AL equations were obtained in [8] using the inverse scattering 
method. In [S] these solutions were derived, using purely algebraic method based on the Fay- 
like identities for the determinants of some special matrices. Here we use notation slightly 
different from one of [9], which makes the following formulae more simple and clear. 



i i 
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The key objects behind the dark-soliton solutions of the ALH are the determinants 

ui {A) = det |I + A\ (4.1) 

with matrices A satisfying 

LA — AR = \ £)(a\. (4.2) 
Here I is the N x N unit matrix, L and R are constant diagonal matrices, 

L = diag (Li, L^) , 

R = dmg(R 1 ,...,R N ), { ' ' 

\i) is a constant iV-column, \l) = {t\, £n) T , and (a\ is a A-row depending on the 
coordinates describing the ALH flows: in our case (a\ = (a (z, z) \ = (a\ (z, z) , a at (z, z)). 
In what follows we use 'shifted' determinants 

uj^ = T^uj, lu^ = T^TjjU; (4-4) 

where 

T 1 ( lj = u(AH l ^j , l = ±l (4.5) 

with 

H^iL-c^iR-ar 1 - (4.6) 

An important property of these determinants, that we repeatedly use below, is the Fay's 
identity 

(f - V) U£ V + (V-C) ^ ^vC + (C - u v = (4.7) 

which can be proved directly. 

Using the limit procedure one can introduce differential operators as 



or 



T^Tf+aw = oj + id 8 c uj + O (8 2 ) (4.8) 

id c A = AX ( (4.9) 

where 

X c = (L - R) (L - CI) -1 (R - (I)' 1 . (4.10) 
One can obtain from ()4.7p many differential Fay's identities of the following type: 

*(C - «)(C - P) D c u a -up = {a-P) [(l^w^) (T c w) - w^] (4.11) 

where 

w a • up = (d^uja) up - co a (d^uip) . (4.12) 
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The matrices L and R used in the ALH context are not independent: they are related 

by 

(L — kT)(R - kT) = —p 2 I. (4.13) 

with constant parameters k and p. Relations of this kind play crucial role in the construction 
of dark solitons for the ALH, so it seems useful introduce the notion of 'duality': two complex 
numbers £ and £* are said to be dual if 

(C-«)(C*-«) = V- (4.14) 

which leads to an alternative definition 

H^H^* = H K , k = const (4.15) 

Many Fay-like formulae can be simplified when rewritten in terms of dual numbers. In 
particular, Eq. (|4.11|) leads to 

iD^ UJ a ■ UJ a * = TjTZ Q *) ( T C 1uJk ) ~ UaUa* (4.16) 

Given some fixed number fi and its dual, which will be denoted by v, v = p*, one can 
construct an doubly infinite set of matrices/determinants 

< = "{K) (4-17) 

where 

A% = AH™H: (4.18) 
and derive from (|4.7p a lot of lattice Fay's identities the most important of which is 

(. ,m\2 _ 2 m—l . m+1 , 2 m m (a i q\ 

with 



Lb — K K — V , 

P» = \~ , Pu = \ • 4.20 

\ p — v \ P — v 

In order to ensure the involution t™ = r~ m (where overline stands for the complex 
conjugation) which appears in physical applications of the Ablowitz-Ladik model one has 
to restrict himself to the case of real k, [i and v, 

v < k < p (4.21) 

(which leads to Im = Im p v = 0) and to choose the matrices L to be of the form 

L = p + y/JT^Dy/jr^H E, E = diag (e*** ) (4.22) 



with real angles ipj (compare with the parametrization (2.27) of the eigenvalues of the 
scattering problem for the ALH used in [8]). 
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Fig. 1. Tp-Tjr* lattice of tau-functions (for n=constant). 



Calculating R from (|4.22|) and (|4.13p one can verify that in this case 



and 



provided 



Am — A-m 
— A n 



A = A. 



(4.23) 
(4.24) 
(4.25) 



The last condition can be met by choosing properly the constants £j and a,- in (|4.2p . By 
straightforward algebra one can get 



Lj - R k = {a-v) 



gi(V>j+i/>fc) 



+ Pll [e^i + e ^*) + 1 



1 + P^e^ 

and rewrite the matrix A, after eliminating excessive constants, as 

A(z,z) jk = D jk c k (z,z) 

where 



D 



( ■•( >S | ; I + Pn COS ' 



(4.26) 



(4.27) 



(4.28) 



and c k (z,z) are some real functions. 



4.2. Dark soliton solutions of ft3.27\ )- {373^) . 

After we have established the relation of our model with the ALH and knowing the structure 
of the ALH dark solitons, we can reformulate the ansatz we use as follows: all tau-functions 
are related by the T-shifts (T u , and T^* for some given £). The sequence p — > r — > a is 
generated by T M , the sequence f — > r — > f is generated by T^*, as is depicted in Fig. 1, while 
the nodes n and n + 1 are related by T„. Thus we can say that our tau-functions occupy 
sites of a three-dimensional lattice. However in what follows we do not use the three-indices 
and adhere to the t™ notation. 
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To find solutions of ()3.2Tj) (|3.32|) we look for our tau- functions in the form 

ra m 2 „n 2 n ,m n n, ,m (a r>rv\ 

T n = Pu. Pu u v w n I 4 - 29 ) 



with similar formulae for f™ and t™ 



where Co™ and are related to uj™ by means of Tp shifts 

4* w ™ 



-m m 2 n 2 ~m~n-m (a or\\ 

T n = Pfi Pu u v u n (4.30) 

T n = Pfi Pv u v U n (4.31) 
nil 



= T e ^ 5 ^ = Tr. 1 ^ (4.32) 



or 

<D™ = T- 1 ^ 1 , C% = T^™"! 1 • (4-33) 

Taking £, 77, £ in (|4.7p being equal to different triples from {//, z/, k, £, £*} one can conclude 
that restrictions 

&V P ( v\ 2 j-v ( ui)\ 2 j-K 

uj £* — n \ vl £* — v \ uv J q — k 



lead to (^29D - (l3T32"1 ) with 





— 1/ 


r 




e 


-A* 


r 





■y A = i — :: 1 * ■■"<■>< 

■> B = M < 4 - 37 > 

Considering the dependence on z and z it can be shown that to meet (|3.27p and (|3.28|) one 
has to take 

id A = AX, —id A = AX (4.39) 

where 

X = const -X K , X = const -X v (4.40) 



We will write the explicit value of the corresponding constants later, after discussing the 
involution (complex conjugation) and reality requirement. 
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4.3. Complex conjugation and parametrization. 

Recalling the definitions of our tau-functions we can present B n and C n as 



where 



1 fli-? 



uv \ £ — a 



n—-k ,~,m— 1 
n-1 



(4.41) 

(4.42) 
(4.43) 

(4.44) 

where overline stands for complex conjugation. A simple analysis of (|4.42|) . (|4.43|) and (|4.18|) 
leads to the restrictions 



C™ = -uv 



u : 



■V-l 



jl — ^* J \£* — V / 

The involution that ensures reality of <r n and that is consistent with the ALH is 

-nm /~i—m 

n n — ~ K - J n 



lm^ = 0, 



i 



and 



H^H^ — H u 

that should be added to the restriction A™ = A~ m discussed above. 



(4.45) 



(4.46) 



The family of £ and £* that satisfy (j4.14f) . ()4.46p and (|4.45p can be parametrized as 

£ =u + ae i(<P+S) 

= v + oe^* - ^ 

where 

P 



(4.47) 



a = \J p — k — v = — 



(4.48) 



and the angles (ft and are related by 

cos (ft = p u cos 6 (4.49) 
Using another parametrization of the matrices L and R, stemming from ()4.22p and (|4.13p . 



Lj = v + a exp (xj + icftj) 
Rj = v + a exp (— %j + i<f>j] 



where the quantities Xj an d <f>j are defined by 



1 



exp (xj +i0j) = — [Pn + e" 

Pv 



(4.50) 



(4.51) 
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one can obtain 

ff € = -diag ( f&-*it ) , H v = -diag ( ^ ( 4 52 ) 

with 

7 ± = - <^ ± - 2 arg [e» - e *(*-fc±«)l . (4.53) 
4.4. Dark solitons of the Toda-Heisenberg chain. 

Now we have all necessary to write dark soliton solutions of the Toda-Heisenberg chain. 
Their structure is given by (|4.42p and (|4.43[) . The dependence on z and z enters through 
the matrices A and the factor uv, 

id A = AX, -id A = AX (4.54) 



where 



or, explicitly, 



Xo 2 

X = —X K , X = \- x (k-v)X v (4.55) 



A ~ A 



-i 



X = (L — R) , X = (L-R). (4.56) 

fi — v (jl — v 

This dependence of A, as follows from the differential Fay's identity (|4.1ip . leads to 

. A (a* — a) r . ,, 

iD oo a * ■ uj a = \lj lo k — u a *u) a \ (4.57) 

ji — v 



and hence to 



A 1 (a* — a) 

iDui a * ■ uj a = [u^uJu — ui a *u} a \ (4.58) 



■n - m+1 ~m _ ^ (£ ~ r-m - m+1 m m+1] 



and 



• T-i ,-,m— 1 ~m ^ (£ f, ,m— 1 , ,m ,~,m— 1 ,-,ml //i 

iDuj n _ 1 -u n = — — — [u n w„_! - w n _ 1 io n \ (4.61) 

,'n, A ,m+l ,-,m _ ^ (C ~ r-m -m+1 m m+1] fi n\ 
w n-l — _ ~ [ UJ n-l UJ n _ w n W n-1 J • 

The extra terms in the right-hand sides of the above formulae can be eliminated by taking 



uv = const x e 



(4.63) 
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with 



ip = X £ — i z + A" 1 tt*—£± z, (4.64) 

LI — V fx — V 

It is easy to verify that all reality conditions are met provided 

|A| = 1. (4.65) 
Now one can check that tau-functions r™ satisfy Eqs. (|3.27p and (|3.28]1 with 



a = A , a = — A (4.6b) 

il — v uv fi — v uv 

Gathering all constants, one comes to the conclusion that for m = the quantities B n , C n 
solve Eqs. (f3TT2D ~ (l3~T5] l with 

r B = r c = -Ae 2iarg( ^ (4.67) 

and 

f B = f c = -A- 1 e- 2iarg(5 *^ ) (4.68) 

Finally, calculating from (|4.47p and (|4.50p the coefficients that describe the z-, z- 
dependence, 

Lj - Rj = 2p, cos e j j+* (4.69) 
ii — v 



and 

" ~ ^ -2ip v sm6 e**, (4.70) 



[X — v 

one can present B n and C n as follows: 



B n (x, y) = B* exp {-2in6 - i<p(x, y)} *"J? ,y \ (4.71) 

An (x, y) 



and 



C n (x,y) = -B n (x,y). (4.72) 

Here 

£* = — |l-p„exp(^)|, (4.73) 

&n(x,y) = det \S jk - D jk exp {2nxk + a k (x,y) +«7^}| , (4.74) 

with the coefficients Dj k and 7^ being defined in (|4.28p . (|4.53[> . The phase <p and the 
functions a k are given by 

tp(x, y) = Ap u sin 9 (x sin «q + V cos «o) + </> (4-75) 
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with 



«o = + arg A 



(4.76) 



and 



with 



dk(x,y) = 4p M cos6» fc (xsina fc + ycosa k ) + a k 



(o) 



(4.77) 



a k = (f) k + arg A (4.78) 

where and aj^ are arbitrary real constants. 

These formulae, together with the vector- matrix correspondence f)2 . 1 [) . lead to the dark 
soliton solutions of the Toda-Heisenberg chain: 



(4.79) 




5. Conclusion. 

We have studied the (2+l)-dimensional system that was reduced to the integrable Ablowitz- 
Ladik equations. Using this reduction we have derived its soliton solutions. It is clear that 
using this approach one can also derive a wide range of other solutions starting from the 
ones already known for the ALH. Thus Eqs. (jl.lip . (|1.9p possess a set of solutions that are 
typical for the integrable systems (solitons, algebro-geometric solutions etc). At the same 
time it is not clear whether this model is integrable or we deal with another example of 
soliton equation that is not integrable (see, e.g., [IQ]), which can occur in multidimensions, 
contrary (as is presumed) to the (l+l)-dimensional case. However, this very interesting 
question, as well as other related questions (such as, e.g., the Painleve test, the symmetry 
analysis), is out of the scope of this paper and may constitute the subject of the subsequent 
studies. 

Appendix A. 

Rewriting the ansatz (|2.15p . (|2.16p in the form 



id B n = Cn-i^n-i (1 — B n C n ) (B n — B n -i) ,^ ^ 

idC n = ( n Y n (1 — B n C n ) (C n+ i — C n ) 



and 

-id B n = ( n Y n (1 — B n C n ) (B n+ i - B n ) 



id C n — Cn-l^n-l (1 — B n C n ) (C n — C n -i) ^ 



one can calculate dd B n as 



dB B n = id {-idB n ) = id( n Y n (1 - B n C n ) (B n+1 - B n ) 



(A.3) 



August 31, 2011 0:16 WSPC/INSTRUCTION FILE the 



Toda-Heisenberg chain 17 



which leads, after repeated usage of (|A. 1[) to 

AnC-n = {id(n + CnCn) Y n (B n+ \ - B n ) — Cn-lCn *n-l (B n — B n -\) . 

Interchanging the d and d derivatives, dd B n = —id (id B n ), one can obtain 

A n £>n = Cn-lCn Y n (B n+ i — B n ) — (id( n -i + Cn-lCn-l) Yn-1 (B n — B n -l) • 

Similar calculations for ddC n lead to 



C 



4 r° 



(—idCn + CnCn) ^n (Cn+1 — Cn) — CnCn-1 Y n -\ (C n — C n -\) 
CnCn-1 Y n (C n+ i — C n ) + (id( n -i — Cn-lCn-l) ^n-1 (Cn ~ C n _i) 



(A.4) 

(A.5) 

(A.6) 
(A.7) 



Comparing the right-hand sides of the above equations with each other and with the right- 
hand sides of Eqs. (|2.10|) one can conclude that conditions 



Cn±l — Cni 
Cn±l = Cnj 

<9Cn = <9Cn = d( n = d( n = 



and 



CnCn 



(A.8) 
(A.9) 
(A.10) 



(A.ll) 



validate the ansatz ([2IES]), (j2T2TI|) . 
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